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ESTIMATORS  OF  MORTALITY  PROBABILITIES 


N.L.  Johnson* 

University  of  North  Carolina  at  Chapel  Hil 


1.  Introduction 


Kuzma  [1967]  and  Prolette  [1975]  have  compared  by  means  of  quite  extensive 


numerical  calculations,  several  estimators  of  q,  the  probability  of  death  in 


an  interval  (given  alive  at  the  beginning  of  the  interval).  It  is  the  purpose 
of  this  note  to  give  some  approximate  algebraic  relations  among  sene  of  these 
estimators.  Thse  confirm,  in  general  terns,  the  conclusions  reached  from  the 
numerical  investigations.  They  also  provide  quite  good  approximations  (to 
four  decimal  places  when  q < 0.15,  apparently)  to  two  rather  complicated  (in 


appearance)  estimators,  in  terms  of  estimators  of  simpler  form.  This  also 


gives  an  improved  appreciation  of  the  nature  of  the  more  complicated,  estimators 


dotation 


Ue  use  the  following  notation  (as  in  Flandt- Johnson  fl°76] ) 


M = number  alive  at  beginning  of  interval 


d = number  of  deaths  in  the  interval 


c = number  "due  to  be  withdrawn'"  in  tr 


interval 


number  of  deaths  among  these  c individuals 


It  will  be  assumed  that  there  are  no  lost  to  follow  up 


or  new  entrants 


during  the  interval  so  that  w+d 
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The  estimators  to  be  discussed  are: 

o/2)  = (N*  - CM ’ 2 - 2Md)1/2}N-1  (sec  ElvebacV  [1950]) 
where  N'  = N - + |d, 

q^  = d (N  - ^w)  2 (the  ’'actuarial'  estimator), 
qC4)  = 1 - j(2M-c)  ~2[{d,2+4(2N-c) (2M-w-2d)}1^“-d’] 2 


q(5)  = d(N  - ic)'1 


(Chians  [1961])  , 
(see  Elveback  [1958]). 


CD 

(2) 

(3) 

(4) 


. •'151  aM) 

The  estimators  q and  q are  much  simpler  in  form  than  q^“  and  q^  ’ . 

3.  The  Approximations 

Using  the  method  outlined  in  the  Appendix,  the  following  annroxinations 
are  obtained. 


where 


5®  * ♦ u ♦ 

(M^r1 


R = N/N ’ = 


N 

■N  - 


q(4)  = U - j(d'/d)q(5))q(5)  . 


(5) 

(6) 
(7) 


Results  of  sone  numerical  calculations  illustrating  the  accuracy  of  these 
approximations,  using  data  from  Table  2 of  Cbianp  [1961]  are  Presented  in 
Table  1. 

For  larger  q,  (5)  may  be  improved  by  addition  of  ^Q^Cq^fl  + ^C3))-l>3^ 
but  values  of  q (>.20,  say)  for  which  this  is  of  any  importance  are  so  large 


that  validity  of  assumptions  made  in  derivation  of  formulae  become  of  over- 
riding importance  (and  usually  uncertain  validity). 


Table  1. 


Approximations  to  various  estimators  of  q 


Interval 

q<2) 

Approx. 

(5) 

-(4) 

q 

Approx. 

(7) 

;(3) 

-CS) 

0-1 

.2434 

.2429 

.2425 

.2426 

.2417 

. 2436 

1-2 

.1822 

.1821 

.1814 

.1815 

.1210 

.1819 

2-3 

.1033 

.1033 

.1030 

. 1030 

.1028 

.1032 

3-4 

.0860 

.0S60 

.9853 

.0858 

.0856 

.0859 

4-5 

.0643 

. 0643 

.0641 

.0641 

.0641 

.0642 

5-6 

. 0584 

.0584 

.0582 

.0582 

.0582 

.0383 

6-7 

.0438 

.0438 

.0438 

.0438 

.0437 

.0438 

7-8 

.0433 

. 0433 

.0432 

.0432 

.0431 

.0433 

8-9 

. 0337 

.0337 

.0337 

.0337 

. 0337 

.0337 

9-10 

.0466 

. 0466 

.0465 

. 04655* 

.04-64 

.0467 

10-11 

.0439 

.0439 

,04J8 

.0438 

.0437 

.0440 

11-12 

.0513 

.0513 

.0511 

.0511 

.0508 

.0513 

* - exactly 


^ f 21 

Except  fcr  the  approximation  to  qv  in  interval  0-1.  t^.ere  is  never  more 
than  0.0001  difference  between  approximation  and  calculated  value  (to  1 
decimal  places).  In  most  cases  the  values  agree  to  4 decimal  places. 

It  is  interesting  to  note  that  for  t?,ese  data  q^  happens  to  be  quite 
a good  approximator  for  q^. 

4.  Pi sous si on 

. Af4') 

From  (7)  it  is  clear  that  the  approximate  value  for  ov  J is  always  less 
than  q^.  This  is  consistent,  with  numerical  results  obtained  hy  Drolette 

Afii  a.  rjy 

[1975].  She  also  finds  that  q w is  greater  than  qv  To  see  that  this  is 


in  accordance  v/ith  our  formulae,  we  note  that  the  right  hand  side  of  (5)  is 


4 


greater  than  if  and  only  if  ?.  > (1  + From  the 

definition  of  P in  (6),  it  follows  that  P.  > (1  + iq^)'1  and  this  is 

1^(31  A("n 

less  than  (1  + , so  it  is  possible  for  the  approximate  value  of  o^“' 

to  be  less  than  q^,  though  the  least  possible  value  O'hen  w = 0)  is 

q(3)[l  - (q<3>)3{«  - 4q<3>  * 2(&h2  * rS^lV1! 


+ 2(ql  + (qw)  } ] 

^ -(3),,  l,-(3).3, 

> qK  •'{1  - lf(q'-  } . 


•'f7')  /-(Ml 

While  the  labor  of  calailating  q and  q w is  not  great,  the  approxi- 
mations put  forward  here  provide  some  insight  into  the  nature  of  these  estima- 
tors, as  well  as  quite  simple  formulae  for  approximate  calculation. 

In  view  of  the  small  differences  usually  found  among  the  numerical  values 

given  by  the  different  formulae,  and  the  usually  dubious  nature  of  assumptions 

„(">!  ~f41 

involved  in  setting  up  the  likelihood  functions  underlying  oK  J and  q v , 

„ f3l  -rn 

there  are  good  grounds  for  using  simpler  formulae  (like  q or  q ) except 


in  very  special  cases. 


APPENDIX 


For  | t | < 1,  we  have  the  expansion 


(l-2t)1/2  . 1 - t - >t2  • (t3  - |t4  ■ It5 

= 1 * |t  - |-(t  + t2  + t3  + t4  + t5)  * 

f i • h - j*a  - 1)_1 . 

1 - (l-2t)1/2  * itn  * (1  - t)-1!  . 


4 3„5 


This  approximation  is  used  in 


q(2'  = (N'/N)  [1  - {1  - 2G;/N’2)d}1/? 


to  yield  (5). 


5 


Similarly,  we  write  (remembering  d ' = c 


~ (4)  . 

J = 1 


4(2N-c) 


^—7  [4(2N-c)2(1  - ^|1)  + 2d'2 


-4d'(2N-c)(l-i£f  ^ ^ l1^2} 

c 4(2N-c) 


and  expand  the  square  root,  to  yield  (7). 
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